Quasi-relativistic behavior of cold atoms in light fields 
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We study the influence of three laser beams on the center of mass motion of cold atoms with 
internal energy levels in a tripod configuration. We show that similar to electrons in graphene the 
atomic motion can be equivalent to the dynamics of ultra-relativistic two-component Dirac fermions. 
We propose and analyze an experimental setup for observing such a quasi-relativistic motion of 
ultracold atoms. We demonstrate that the atoms can experience negative refraction and focussing 
by Veselago-type lenses. We also show how the chiral nature of the atomic motion manifests itself 
as an oscillation of the atomic internal state population which depends strongly on the direction of 
the center of mass motion. For certain directions an atom remains in its initial state, whereas for 
other directions the populations undergo oscillations between a pair of internal states. 

PACS numbers: 03.75.-b, 39.25.+k, 42.50.Gy, 42.25.Bs 



INTRODUCTION 

Two-dimensional (2D) quantum systems are a source 
of many remarkable phenomena. A striking example in 
this respect is provided by the properties of electrons 
in graphene [11,11, H, 0, H, H, 0, B Q — a two dimen- 
sional hexagonal crystal of carbon atoms. Near the Fermi 
level the electrons in graphene behave like massless ultra- 
relativistic two-component Dirac fermions [lo| moving 
with a velocity vp which does not depend on the mo- 
mentum. This leads to a number of distinct effects in 
graphene and graphene bilayers, such as a half integer 
quantum Hall effect 0, 0, Il| and the Klein paradox [sj . 
Furthermore it has been suggested recently that electrons 
in graphene should exhibit a negative refraction [^, Q at 
a potential barrier, similar to the electromagnetic waves 
refracting in a counterintuitive way at an interface with 
a material characterized by a negative dielectric permi- 



tivity and a negative magnetic permeability [11 



In this Rapid Communication we show how cold atoms 
obtain ultra-relativistic properties of Dirac fermions if 
manipulated by laser beams. We suggest to use three 
laser beams acting on atoms in a tripod configuration 

. The light beams induce an effective vector potential 
(the Mead-Berry connection 1^, which influences 

the atomic center of mass motion [ij]. We demonstrate 
that by choosing proper light fields the vector potential 
can be made proportional to an operator of spin 1/2. For 
small momenta the atomic motion becomes equivalent 
to the ultra-relativistic motion of two-component Dirac 
fermions, as is the case for electrons in graphene. We 
propose and analyze an experimental setup for observing 
such a quasi-relativistic behavior of the cold atoms. We 
show that the atoms can experience negative refraction 
and focussing by Veselago-type lenses |11|. 
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FIG. 1: (Color online) An atom interacting with three laser 
fields in a tripod configuration. 



Interestingly the velocity of the quasi-relativistic atoms 
is of the order of a centimeter per second. This is ten 
orders of magnitude smaller than the speed of light in 
vacuum c ~ 3 x lO^m/s. For comparison, the velocity 
of the Dirac-type electrons in graphene vp ~ 10^ m/s is 
only three hundred times smaller than c Thus the 
ultra-relativistic behavior of cold atoms manifests itself 
at extremely small velocities. Note also that our proposal 
does not need a lattice and thus operates in a continu- 
ous regime in contrast to recent papers on the quasi- 
relativistic dynamics of cold atoms in one-dimensional 
[ist or 2D hexagonal (graphene-type) fl^ lattices. 



FORMULATION 

Let us consider the adiabatic motion of atoms in the 
field of three stationary laser beams. The beams couple 



four internal atomic levels in a tripod configuration |13l . 
14l | , in which the atoms are characterized by three lower 
levels |1), 1 2) and |3) and an excited level |0) shown in 
Fig. [1^ [20[. The j-th laser induces a transition (with 
a Rabi frequency ^j) between the j-th lowest level and 
the excited level |0). We shall be interested in a scheme 
where the first two lasers have the same intensities and 
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counter-propagate in the x-direction while the third one 
propagates in the negative y-direction (see Fig. [TJs), i.e. 
rii = nsmee-"^'^/V^, = »sin6>e^''"'/V2 and f^g = 
ncosee-"'y, where n = ^|f7i|2 + If^aP + jr^gP is the 
total Rabi frequency, and the mixing angle defines the 
relative intensity. 

The electronic Hamiltonian of the tripod atom reads 
in the interaction representation 14 1 



i?o = -n(ni\0){l\ + r!2|0)(2| + r!3|0)(3|j + H.c. . (1) 

The Hamiltonian Hq has two eigen-states of zero energy 
containing no contribution from the excited state |0): 
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The states |Z?i) and \D2) known as the dark states 
141 depend on the atomic position through the spatial 
dependence of the Rabi- frequencies Vlj . 

The adiabatic approximation is carried out neglect- 
ing transitions from the dark states to the bright state 
\B) -^Vll\\)+Vll\2)+D.l\i). The latter is coupled to the 
excited atomic state |0) with the Rabi frequency 17, so the 
two states \B) and |0) split into a doublet separated from 
the dark states by the energies ztil. The adiabatic ap- 
proximation is justified if f2 is sufficiently large compared 
to the two-photon detuning due to the laser mismatch 
and/or Doppler shift. In that case the internal state of 
an atom evolves within the dark state manifold. The 
atomic state- vector |<i>) can then be expanded in terms 
of the dark states according to 1$) = (r)), 
where vE'j (r) is the wave-function for the center of mass 
motion of the atom in the j-th dark state. 

Thus atomic center of mass motion is described by 
a two-component wave-function — (5'i,^2)^- The 
column-matrix 5" obeys the Schrodinger equation [l^ 
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-ihV ~Af + V + ^ 



(4) 



where m is the atomic mass, and A, <i> and V are 2x2 
matrices. The gauge potentials A and $ emerge due 
to the spatial dependence of the atomic dark states. 
The reduced 2x2 matrix A with the elements A„ m = 
ift(_D„(r)|V£'m(r)) represents the effective vector poten- 
tial known as the Mead-Berry connection 14, 15, 3) 12 1- 
The 2x2 matrix <I> acts as an effective scalar potential. 
The external potential V confines the motion of the dark- 
state atoms to a finite region in space. Specifically we 
have V^,n = {Dn{r)\V\D^{r)) with V = yi(r)|l)(l| + 
V2{v)\2){2\ + y2(r)|3)(3|, where V,(r) is the trapping po- 
tential for an atom in the internal state j — 1,2,3. Note 
that the potential Vj can also accommodate a possible de- 
tuning of the j'-th laser from the resonance of the j — s- 
transition. 



ATOMIC MOTION IN THE FIELD OF THREE 
PLANE WAVES 

The potentials A, <I> and V have been considered in 



Ref. [1J| for arbitrary light fields. In the present config- 



uration of the light fields, the potentials take the form 
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(5) 
, (6) 
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where the external trapping potential is assumed to be 
the same for the first two atomic states, V\ = Vi- 

In what follows we take Vi, — V\ = h^n^ s\v?(0)/2m. 
This can be achieved by detuning the third laser from 
the two-photon resonance by the frequency Awa = 
Kk^ SIT? 6 /2m. Thus the overall trapping potential sim- 
plifies to V + ^ = Vil (up to a constant), where I is the 
unit matrix. In other words, both dark states are affected 
by the same trapping potential Vi = Vi (r) . 

Furthermore we take the mixing angle 9 = 6^) io be 
such that sin^6'o = 2cos9q, giving cos^o = \/2 ~ 1- In 
that case the vector potential can be represented in a 
symmetric way in terms of the Pauli matrices Gx and cr^ 



A = ftK'(-e^o-^ + eydz) + hK^eyl, 



(8) 



where k' = k cos 6*0 ~ 0.414k and ko = '*(! ~ cos 6*0). 
Although the vector potential is constant, it can not 
be eliminated via a gauge transformation, because the 
Cartesian components Ax and Ay do not commute. Thus 
the light-induced potential A is non-Abelian. Note that 
non-Abelian gauge potentials can also be induced in op- 
tical lattices using other techniques 21 1. 

It is convenient to introduce new dark states: 
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\D[) = _(|i?l)+^|i?2))e-«^ 



\D'\ 
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{\D{)-i\D2))e''''y. 



(9) 
(10) 



The transformed two component wave-function is 
related to the original one according to — 
exp(— iKoy) exp (— i^CTa;) 5". The exponential factor 
exp(— iKo?/) induces a shift in the origin of the momentum 
k —t k — KflGy. With the new set of dark states we get the 
vector potential A' = —hn'tj^, where a±^ — Bx^^x + ^y^y 
is the spin 1/2 operator in the xy plane. The transformed 
equation of the atomic motion takes the form 
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(11) 



In this way the vector potential governing the atomic 
motion is proportional to the spin operator <t± . 
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If the trapping potential Vi is constant, we can consider 
plane- wave solutions, 



*2k 



(12) 



where is the eigen-frequency. The k-dependent spinor 
obeys the stationary Schrodinger equation i/k^k = 
^t^k^k, with the following k-dependent Hamiltonian 
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(k + K'a^^f + Vi. 



(13) 



For small wave- vectors (fc <C k') the atomic Hamiltonian 
reduces to the Hamiltonian for the 2D relativistic motion 
of a two-component massless particle of the Dirac type 
known also as the Weyl equation (2^ . 



-ffk = fivok ■ <T± +Vi 



(14) 



where vq = hn'/m is the velocity of such a quasi- 
relativistic particle. The velocity vq represents the re- 
coil velocity corresponding to the wave- vector a' , and is 
typically of the order of a centimeter per second. 

The Hamiltonian iJk commutes with the 2D chirality 
operator dk = k • cr^/fc. The later has the eigenstates 
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(15) 



for which crk^'k = ^l^^k ■ Here the chirality is associated 
with the subspace of the dark states defined by the basis 
in Eqs. ([9|l~ (fT0|) rather than with the spin states in the 
usual sense. The corresponding dynamics can therefore 
also be illustrated by for instance a Poincare sphere. 

The eigenstates (ITSt are also eigenstates of the Hamil- 
tonian i?k with eigen-frequencies Wk = w^. In what 
follows the atomic motion is assumed to be confined in 
the xy plane. The dispersion then reads 



hvo{k^/2K' ±k) + Vi + 



mvl , 



(16) 



see Fig. [2] in which Vi — —mv'^12. The atomic mo- 
tion in different dispersion branches is characterized by 
opposite chirality if the direction k/fc is fixed. For 
small wave-vectors (fc ^ k') the dispersion simplifies to 
tvjj^ — zthvok + Vi + toUq, where the upper (lower) sign 
corresponds to a linear cone with a positive (negative) 
group velocity, = ±vo- Exactly the same dispersion 
is featured for electrons near the Fermi level in graphene 

UMMMM- 



PROPOSED EXPERIMENT 

To observe the quasi-relativistic behavior of cold 
atoms, we propose the following experimental situation. 
We suppose initially an atom (or a dilute atomic cloud) 




FIG. 2; (Color online) Upper (red solid) and lower (green 
dashed) dispersion branches for a tripod atom in light fields 



is in the internal state |3) with a translational motion 
described by a wave-packet with a central wave-vector 
kin and a wave- vector spread Ak ^ fcjn, i.e. I^'in) = 
'0(r)e''''° '"|3), where the envelope function tp{r) varies 
slowly within the wavelength Ain = 2?:/ kin- The cold 
atoms can be set in motion using various techniques, e.g. 
by means of the two-photon scattering which induces a 
recoil momentum Kk-in = h]s.2phot to the atoms, where 
k2phot is a wave- vector of the two-photon mismatch [23| . 

Initially all three lasers are off. Subsequently the lasers 
are switched on in a counterintuitive manner, switching 
the lasers 1 and 2 on first followed by the laser 3. At the 
beginning of this stage the internal state |3) coincides 
with the dark state I-D2), so the original and transformed 
multi-component wave-functions are given by: 



V2 \ 1 



(17) 

where the transformed wave-function ^' is obtained ex- 
pressing the original dark state ID2) through the new 
ones \D'i 2), with k ~ kjn — KpSy being a new central 
wave-vector. If the laser 3 is switched on sufficiently 
slowly, the atom remains in the dark state I-D2) dur- 
ing the whole switch-on stage. Yet the duration of the 
switching on should be short enough to prevent the dy- 
namics of the atomic center of mass at this stage. To 
have ultra-relativistic behaving atoms, the wave-number 
k should be small A: <C k, so that k is a small contribu- 
tion to kin = KQey + k. In addition, the wave-number 
spread Ak k, i.e. the width of the atomic wave-packet 
is much larger than the central wave-length. The subse- 
quent centre of mass motion of atoms in the laser fields 
is sensitive to the direction of the wave-vector k. 
(i) If k = zLkey, the wave-function ([T7)) reads: 
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(18) 
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FIG. 3: (Color online) Negative refraction of cold atoms at 
a potential barrier. The incoming and outgoing atoms are in 
the upper (red) dispersion branch, whereas the atoms inside 
the barrier are in the lower (green) one. 



The upper (lower) sign in k = ±key corresponds to a sit- 
uation where the atom is characterized by a positive (neg- 
ative) chirality, hence being in the upper (lower) disper- 
sion branch. In both cases the atomic wave-packet prop- 
agates along the y axis with the velocity vq = Byhn' /m. 

(ii) If the wave- vector is along the x axis (k — ke^), 
the multi-component wave-function (|17p takes the form 
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1 

±1 



(19) 

where c± = {—i ± l)/2. In that case the initial wave- 
packet splits into two with equal weights (|cj_ | = 1/2) and 
the same wave-vector k. The wave-packets are charac- 
terized by different chiralities and thus move in opposite 
directions with the velocities vq = dze^hn' /m. 

Suppose the time is sufficiently small {vot < d), so the 
wave-packets of the width d are not yet spatially sepa- 
rated. The internal atomic state will then undergo tem- 
poral oscillations between the dark states I-D2) and \Di) 
with a frequency equal to ui^ — uj^ = 2vok. Such an in- 
ternal dynamics can be detected by switching the laser 3 
off at a certain time. This transforms the dark state ID2) 
to the physical state |3). Subsequently one can measure 
population of the state |3) for various delay times and 
various wave-vectors k. The chiral nature of the atomic 
motion will manifest itself in the oscillations of the popu- 
lation of the atomic state |3) if k is along the x axis, and 
the absence of such oscillations if k is along the y axis. 

Furthermore as a consequence of the constructed 
Hamiltonian (|14p . the quasi-relativistic atoms can show 
negative refraction at a potential barrier and thus exhibit 
focussing by Veselago-type lenses [111, Il2| . Consider inci- 
dent atoms that are in the upper dispersion branch and 
propagate along the y axis with a wave- vector k = key. 
Let us place a potential barrier of a height 2hvok at an 
angle of incidence a (see Fig. |3]). Inside the barrier the 
atoms are transferred to the lower dispersion branch with 
kj — — A;[cos(2Q!)ej^ -I- sin(2Q;)ea,]. This would lead to the 



negative refraction of cold atoms at the barrier as shown 
in Fig. [3] Thus the potential barrier can act as a flat lens 
which refocuses the atomic wave-packet. 

In summary we have shown how the atomic motion 
can be equivalent to the dynamics of ultra-relativistic 
(massless) two-component Dirac fermions. As a result 
the ultracold atoms can experience negative refraction 
and focusing by Veselago-type lenses. In addition, we 
have investigated another manifestation of the chiral na- 
ture of the atomic motion through dynamics of the pop- 
ulation of the internal atomic states. 
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